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Abstract 

Recent supernovae (SNe) detections have motivated renewed interest in SN shock breakouts from 
stars surrounded by thick winds, including predictions of observable hard X-rays. Wind breakouts on 
timescales of a day or longer are currently the most probable for detection. Here we study the signal 
that follows such events, assuming a wind density profile oc r~^, starting from the breakout of the 
radiation mediated shock and tracing the evolution of the collisionless shock which forms afterwards. 
The emission contains two spectral components - soft (optical/UV) and hard (X-rays and possibly soft 
7-rays). We find that during the breakout, the soft component temperature can vary significantly from 
one event to another (lO"* — 10^ K), where events with longer breakout time, iho, are generally softer. 
The hard component is always a minute fraction, ~ 10""*, of the breakout emission, and its fraction 
of the total luminosity rises quickly afterwards, gaining dominance at 10 — 50 ibo- The spectral 
evolution of the soft and hard components, as well as the prospects for detection of X-rays, depend 
mostly on the breakout time. In early breakouts {tho ^ 20 d for typical parameters) both components 
become harder after the breakout, and the hard component becomes dominant while the luminosity is 
still comparable to the breakout luminosity. In late breakouts {tho ^ 80 d for typical parameters) the 
soft component becomes softer with time and the hard component becomes dominant only after the 
luminosity has dropped significantly. In terms of prospects for X-ray and soft gamma-ray detections, 
it is best to observe 100-500 days after explosions with breakout timescales between a week and a 
month. 



1. INTRODUCTION 



The breakout of a supernova (SN) shock through the 
stellar surface has been an active field of anal ytic and 
nume ri cal rese a rch fo r several dec a des (e.g. iColgajd 
1971 IWeaved [l976l iFaia ITOTSL IKlein fc Chevalier 



19781 llmshennik et all 119811 lEnsman fc BlITrows ,1992. 
Matzner fc McKed Il999f ) Recent advancement of 



observational facilities led to the discovery of several 



20061 iSoderberg et al. 


120081 IGezari et al.l 20081: 


Schawinski et al. 20081: 


Modiaz et al.l 20091 Ofek et al.l 


20101 IGezari et al.l 120101 


Arcavi et al.l 120 111) which 



motiv ated a rev isiting of the theor y of shock br e akout 
fe.g.. iChevalie r & Franssoi ] [2001 iKatz et al.l [20Tol: 



Piro et al 
20Tla[ 



201ir 



2010; Nakar & Sari' '2010; 'Katz et al.' 

Rabinak & Waxman 2011; Couch et al. 

.Nakar fc Saril l20Tlt [Ch evalier fc Irwiij l20TlI 

Balberg fc Loebll2011l : lMoriva fc Tominagall201lD 

The main observational challenge posed by shock 
breakouts from stellar surfaces is their short duration (no 
more than an hour even in the case of a red supergiant). 
However, in some cases massive mass loss prior to the 
SN explosion can extend the bre akout signal, facilitating 
its de tection. lOfek et all (|2010D and IChevalier fc IrwinI 
(|2011h show that a breakout through a thick wind may 
occasionally extend long enough to account for the com- 
plete SN luminosity light curve, and suggest that at least 
some type Iln SNe are in fact powered by such breakouts. 

When a star goes through a SN explosion, a radiation 
mediated shock (RMS) is driven through its envelope, ac- 
celerating through the sharp density drop near the stel- 



lar edge (jSakurail I1960D . If the star is not surrounded 
by a thick wind then the shock breaks out of the stellar 
surface, producing an intense short pulse of UV- X-ray 
photons. If, on the other hand, the star is surrounded 
by a wind with an optical depth between the stellar edge 
and the observer > c/v, where c is the speed of light and 
V is the shock speed, then the shock continues into the 
wind without releasing photons to the observer. At this 
point the shock starts decelerating and a reverse-forward 
shock structure is formed. The reverse shock is driven 
into the SN ejecta, which is characterized by a sharp rise 
in the amount of energy that is carried by slower mov- 
ing material. The reverse shock is driving an increasing 
amount of energy into the shocked region, which keeps 
accumulating until the wind optical depth drops to c/v. 
Beyond this point the wind cannot contain anymore the 
radiation mediated forward shock (the shock width be- 
comes ^ R), and the shock is breaking out of the wind, 
releasing all its energy as an intense pulse. This pulse 
can contain much more energy than a wind-less break- 
out, over timescales that are much longer. The general 
pr operties of a wind breakout were recentl y disc ussed 
bv lOfek et all (|20Toh and IChevalier fc IrwinI jMp), who 
find the pulse luminosity and tem perature assum ing ther- 
mal equilibrium, and by Balberg fc Loebl ()201lF ) who also 
show how the progenitor mass and explosion energy can 
be constrained by the pulse observational properties. 

If the wind density does not fall abruptly, the energy 
behind the shock keeps growing after breakout, and the 
characteristic evolution of the resultant signal supplies 
an important observational probe for the interaction of 
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the forward shock with the wind. iKatz et all ()2011bO 
showed that if the breakout of the RMS occurs within 
the wind, i.e. the matter consisting the breakout layer 
is spread over a scale comparable to the breakout ra- 
dius rather than being concentrated within a thin layer, 
then a collisionless shock is bound to develop and re- 
place the RMS in accelerating the matter ahead. This 
collisionless shock heats the immediate pos t-shock elec- 
trons to Te ~ 60 keV (jKatz et al.ll2011b[) . where Te is 
the electron temperature and ks is the Boltzmann con- 
stant. The cooling processes of these electrons affect the 
observed temperature (and possibly the luminosity evo- 
lution), which is further influenced by the interplay be- 
tween photons and electrons along the photons diffusion 
path to the observer. 

Here we examine the evolution of the luminosity and 
the observed temperature at and following a wind shock 
breakout. We consider only breakouts through a wind 
that is dense enough to produce a breakout pulse lasting 
days or longer, which are easier to detect. We also re- 
strict the discussion here to shocks that are cooling fast 
during the entire evolution, i.e., the shock heated plasma 
cool within a dynamical time. We show that this con- 
dition is satisfied for regular SNe parameters. We pay 
attention to deviation from thermal equilibrium and to 
hydrodynamic, diffusion and cooling timescales involved, 
and show that the breakout timescale is the most domi- 
nant parameter in determining the post-breakout evolu- 
tion. We map the different types of post breakout behav- 
iors and describe their light curves an d spectra evolution. 

As t his work was near completion, iChevalier fc IrwinI 
()2012l) submitted a paper, presenting some of the conclu- 
sions derived below in the context of explaining the low 
X-ray luminosity of SN 2006gy, observed by Chandra at 
an age of 3-4 months. The low luminosity is explained 
as the result of three factors: two, which are discussed 
in details below, are the domination of inverse Compton 
over free-free emission in cooling the hot post-shock elec- 
trons, and the energy degradation of the hard free-free 
emitted photons over the diffusion path. The third is 
photoabsorption of soft X-rays by partially ionized pre- 
shock gas. This latter process is not considered in this 
work. 

2. HYDRODYNAMIC EVOLUTION AND 
CHARACTERISTIC SCALES 

We consider the interaction of SN ejecta with a dense 
stellar wind that follows a standard profile pw = Dr~^, 
where D is a constant and r is the radius. This wind pro- 
file implies r cx r~^, where r is the optical depth of the 
wind at the location of the shock. It also implies a photon 
diffusion time scale to the observer, that is independent 
of the radius (up to a logarithmic factor), tdiff ~ kD/c, 
where k is the cross section per unit of mass (throughout 
the paper we assume that the gas which dominates the 
optical depth is fully ionized, and take k = 0.34 cm^g~^). 

During the interaction, a forward shock is driven into 
the wind while a reverse shock is driven into the ejecta. 
iChevalieil ()1982fl derives a self similar solution of the 
forward-reverse shock structure, assuming that both the 
ejecta and the wind density profiles are power-laws: 

r{t) cx (1) 
for pu, cx r^^ and pej cx v~". The ejecta density profile 



through which the reverse shock is propagating is deter- 
mined by the SN shock that crosses the stellar envelope. 
The density profile of the fastest moving ejecta can be 
approximated by a power-law with an indeju n ^ 10 — 12 
for all progenitor types. The ejecta density profile flat- 
tens towards the slower layers which contain most of the 
ejecta mass and energy, (e.g. iMatzner fc McKee| [l999). 
As the reverse shock approaches these layers the evolu- 
tion is no longer strictly self-similar, but it can be ap- 
proximated by the self-similar solution using n — 7 (e.g., 
IChcvalicr & Irwin 2011, 2012). Given that during the 
reverse-forward shock interaction phase n > 7 we here- 
after approximate « 1 during this phase. Finally, 
once the reverse shock ends crossing the ejecta only the 
forward shock remains, entering a Sedov- Taylor phase if 
the shock is adiabatic or a snowplow phase if it is radia- 
tive. 

The density power-law indices of the ejecta and the 
wind are set by the pre-explosion stellar evolution and 
are not expected to vary much between various SN pro- 
genitors. The difference in the hydrodynamical evo- 
lution during the forward-reverse shock phase between 
various SNe is set by the two normalization factors 
of these density profiles. These could be determined 
by two observables such as the duration of the break- 
out pulse, th n, which is set by the wind density only 
thn ^ /tZJ/clOfek et all IChev alier fc Irwid [20Tll: 

iBalberg fc Loebll20lH ). and the breakout luminosity Lto- 
Below we present the evolution during the interaction 
phase using tto and the shock breakout velocity Vbo , (in- 
stead of Lbo) which is related to the observables via 

Vbo = {-^-^Lbo^^f^ (IBalberg fc Loebl[20ll . The evo- 
lution after the interaction phase ends depends on three 
parameters and we present it using tbo, the total explo- 
sion energy, E, and ejecta mass, M^j. 

The breakout from the wind takes place once Tbo ~ 
c/vbo, where r^o is the wind optical depth to the observer 
at the breakout radius. This sets the mass swept by the 
shock at the time of the breakout: 
47rc 

rubo ~ Vbotla w 5 X 10^^ MQVbo.gtlo d > (2) 

K 

where Vbo,9 = loa^'cm/s ^^'^ ^bo,d = Yd^- During the 
ejecta-wind interaction the wind mass collected by the 
forward shock is comparable to the ejecta mass collected 
by the reverse shock. Therefore, using the velocity pro- 
file of an ejecta re leased by a SN explosion (e.g., from 
INakar & Sari [201(11 : hereafter NSIO), and requiring that 
a wind mass collected by the forward shock until the 
breakout, nibo, is accelerated to Vbo by the explosion, im- 
plies: 

« 10« cm/s M,-r'Elft^^, (3) 

where A/io — and Er,i — jqst-^- This equation is 

applicable until mbo = M^j. It implies that for typical 
SNe (e.g., not extremely high or low E^i/Mio) and in 

_ n + l 

The mass profile for m <C M^j is m oc ?; o.i9n wfiere n (not 
to be confused witii n) is tfie power-law index of the pre-explosion 
stellar density profile near the edge (1.5 for a fully convective enve- 
lope and 3 for a fully radiative one). This profile translates, under 

free expansion, to pgj oc t '^v o.Wn , or n 10 for a radiative 
envelope and n 12 for a convective one 
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the range of tbo that we consider (> 1 day) the break- 
out velocity does not vary much. Thus, we consider here 
only a narrow range of Vbo,9 = 0.3 — 1.5. The upper 
limit of this range ensures that shocks that follow break- 
outs with timescales of days or longer remain fast cooling 
(i.e., the shock heated plasma radiate its energy within 
a dynamical time) at all times (see section 13 3p . The 
lower velocity limit allows us to neglect recombination 
effects during the breakout. Interestingly, Equation [3] is 
independent of the progenitor radius and is very weekly 
sensitive to the exact density profile near the stellar edge. 
It is therefore applicable to red and blue supergiant as 
well as Wolf-Rayet progenitors. Thus, if the shock is 
breaking out from a thick wind, the wind density is the 
parameter that affects the observed radiation the most 
and the progenitor type cannot be identified based on 
the light curve alone. 

An important time scale, which affects the hydrody- 
namic evolution, is the time at which the wind mass col- 
lected by the forward shock, m, is comparable to the total 
ejecta mass. At that point the reverse shock stops play- 
ing a role and the forward shock is slowing down rapidly. 
The deceleration rate depends on whether the shock is 
adiabatic, or a significant fraction of the post-shock in- 
ternal energy is radiated away. If it is adiabatic then the 
evolution follows the Sedov- Taylor self similar solution. 
If, on the other hand, all the internal energy behind the 
shock is radiated within a dynamical time scale, then 
the shock enters the so called snowplow phase. We con- 
sider only fast cooling shocks, which are in the snowplow 
regime when m > Mej and t > t^o- 

Next we find the time of transition to snowplow, tsp, 
in case that tbo < tsp- Equation [T] for r{t) can also be 
applied for the snowplow phase using n — 4 (the Sedov- 
Taylor phase matches n = 5)0. Therefore, in case that 
tbo <tsp the swept mass after breakout accumulates as 



takes place before tsp then 



m{t) 



rribo-, 



thn<t <t 



SP 



(tsp) 



1/2 



tsp < t 



(4) 



Since the snowplow regime begins once m = A/gj, it fol- 



lows from Equation [2] that if 4o < SOAfiV^^'si d then 
tbo < tsPj where we approximate the shock velocity at 
the transition to the snowplow phase with ^JIE jMp^y 
The transition time to snowplow phase occurs then at 



tsp 



7x10^ CXoX'^'d 



(5) 



If m = Mgj before breakout, then the evolution follows 
a Sedov- Taylor expansion prior to breakout and assumes 
a snowplow evolution afterwards. In such case we say 
that tsp < tbo and the mass behind the forward shock 
accumulates as 



m(t) ss nibo 



t 

tbo 



1/2 



tsP < tbo < t 



(6) 



3. TEMPERATURE AND LUMINOSITY EVOLUTION 
The bolometric luminosity of a fast cooling shock is 
L « 2TTr'^pv^ (X (X t~"^ . Therefore, if the breakout 

^ Note that using n = 4 in Equation [T]provides the cor rect evo- 
lution in the snowplow regime, although IChevalierl 119821 ) derived 
it for the interaction phase only. 



Lit) 



5 X 10 



43 



tbo.dvl^g (t^) 



-0.3 



^ t « tsP 

1.5 



2 X m^Hbo,,M-,'/'ElC ■' ^ tsp < t 

(7) 

while if tbo > tsp then 

L{t) - 5 X 10^3 i^^^^ y3^fl_\ ^ . tgp < < t. 

' \tbo J S 

(8) 

Note that once the radiation escape time becomes shorter 
than the dynamical time (i.e., after breakout), then as 
long as the shock is cooling fast, the bolometric luminos- 
ity is simply the kinetic energy flux through the shock. 
In contrast to regular SNe, it is independent of the gas 
opacity, and thus insensitive to processes such as recom- 
bination or to the gas metallicity. 

Following breakout, the collisionless shock heats pro- 
tons. The protons then transfer heat to electrons, 
which are cooled by radiation. The balance between 
the electrons heating and cooling rates sets the electron 
temperatur e, Tp, d uring the post-shock cooling process 
pCatz et ani2011bl ). The breakout balance temperature 
is ^ 60 keV, independent of the shock velocity. If the pro- 
tons temperature, kBTp k, j^rUpV^ k, 200wg keV where 
mp is the proton mass, is lower than the implied balance 
temperature, then T,^ ^ Tp. Thu^ 



B-Le,bo 



200v, 



ho. 9' 



60] keV. 



(9) 



It is straight forward to show that after breakout Te ap- 
proaches Tp (if they where not similar already at break- 
out). Now, since we consider Vbo,9 ~ 0.3 — 1.5, the elec- 
tron and proton temperatures are not very different al- 
ready during the breakout, and at i < tsp the proton 
and electron temperatures are approximately constant. 
At t > tsp the proton temperature drops, and Tg follows 
it. Therefore, 



ksTeit) 



60 keV 

20E,,M^,' (^) 



keV 



t < tsP 
t > tsP 



(10) 

While the electrons behind the collisionless shock are 
hot, the observed radiat ion at the breakout, aii d soon 
after it, is much softer (jChevalier fc IrwinI 120121) . The 
reason is twofold. First, at breakout and for some time 
after it, the main cooling source is inverse Compton (IC) 
over optical-UV photons. These photons are generated 
by the unshocked wind that is heated by the outgoing 
radiation. Thus, the energy in hard photons, which are 
emitted by free-free emission, is suppressed by the free- 
free to IC cooling rate ratio. Second, even the small num- 
ber of emitted hard photons do not reach the observer 
directly. Their energy is reduced via Compton scattering 
on the unshocked, lower temperature, wind electrons. 

^ This equation assumes that the shock does not couple the elec- 
trons to the protons and that the electron heating is by Coulomb 
interactions. This is not necessarily true, in which case electrons 
are heated by the shock to Tp and then cool down until the radia- 
tive cooling and Coulomb heating is balanced. Since the difference 
between Tp and Te is not large under any of the two scenarios their 
predictions for the observed signature are similar. 
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Thus, the observed spectrum will have two compo- 
nents, a soft one and a hard one. The softer component, 
which dominates the luminosity at early times, is pro- 
duced by IC and its temperature is determined, mostly, 
by the deviation of the photon field in the unshocked 
wind from thermal equilibrium. The harder component 
is emitted by free-free. Its temperature is determined by 
the energy loss of a single photon while it diffuses out of 
the wind and its fraction of the total luminosity is set 
by the product of this energy loss and the free-free to IC 
cooling rate ratio. Below we calculate the evolution of 
the temperature of the soft component, Tso/t, as well as 
the temperature and luminosity, Thard and Lhard- 

3.1. Soft component 

The shock-heated electrons and protons are confined 
to a narrow layer behind the shock. Although they are 
the source of the observed luminosity, they do not neces- 
sarily set the typical photon temperature. This temper- 
ature is set during the diffusion of the photons towards 
the observer through the unshocked wind, as long as it 
is optically thicl|3. The unshocked wind is heated by the 
diffusing radiation, which produces photons, mostly by 
free- free, which share energy with the diffusing radiation. 
These photons also diffuse back into the hot layer behind 
the shock, serving as a cooling source. The soft compo- 
nent temperature is therefore set by the ability of the 
unshocked wind to produce photons over the radiation 
diffusion time, tdi// ^ tbo- The radiation energy den- 
sity at any point along the diffusion path is trad ~ 

We define Tbb = {^rad/cY^^ , where a is the radiation 
constant. If photons are abundant enough and the ra- 
diation is in thermal equilibrium then the photon tem- 
perature at the shock is Tbb- The observer is seeing, 
however, a lower temperature. The reason is that ther- 
mal equilibrium is also kept at radii that are larger than 
the shock radius, where the radiation energy density is 
lower. Thus, the observed temperature is set by the radi- 
ation energy density at the largest radius at which ther- 
mal equilibrium is kept. If, in contrast, the photons at 
the shock are out of thermal equilibrium, then the ob- 
served temperature can be much higher than Tbb, and 
is set by the number of photons that are generated over 
the diffusion time. This process is the same one as in the 
case of a shock breakout from a stellar surface where no 
wind is present, a nd is discussed in length in NSIO and 
iKatz et al.l (|2010l ). 

Assuming that photons and electrons in the unshocked 
wind have the same temperature, the observed temper- 
ature can be determined using the coupling coefficient 
r/ EE 5 X lQ^^^T^lp~H~ljf defined in NSIO, where aU 
quantities are in c.g.s. units. Here rj is the ratio of the 

Fast cooling dictates that the shocked plasma is cooling, form- 
ing a dense cold layer behind the hot shock heated plasma. This 
cold layer does not affect the observed soft component. The rea- 
son is that after the breakout the photons escape on a time that 
is shorter than the dynamical time, and therefore the energy den- 
sity behind the shock is not dominated by radiation. Since the 
shock is radiative the shocked plasma is compressed to keep the 
pressure balance and once free-free cooling becomes dominant, its 
cooling rate increases with the compression and a run-away cool- 
ing decouples the plasma temperature from the rest of the system, 
preventing it from contributing a significant number of photons to 
the radiation field that cools the shock heated plasma. 




0.3 0.5 0.7 0.9 1.1 1.3 1.5 

V [10^cm/s] 



Fig. 1. — The observed breakout temperature as a function of the 
breakout velocity for 1^,^ = 20 d, demonstrating the strong effect of 
deviation from thermal equilibrium on the observed temperature. 
The relatively small range of relevant velocities 0.3 < i>i,o,9 < 1-5 
corresponds to a range spanning over 2.5 orders of magnitude in 
temperature, lO^K < Tsoft,bo < 5 X 10^'K. 

photon density needed to maintain thermal equilibrium 
to the photons produced per unit of volume during the 
diffusion time, assuming that the electrons' temperature 
is Tbb- When its value at the shock radius is rys < 1, 
the radiation is in thermal equilibrium and the temper- 
ature is set by the electrons at the radius where = 10. 
When > 1, the observed temperature is higher than 
Tbb and is calculated as described in NSIO (with some 
modifications described below). 
At breakout, aT^j^/3 = (6/7)p bovlg, implying 

and yielding 

ribo ~ ^vLXi^- (12) 

At any given time rj cx r^^/^, implying that if 77^0 ^ 1, 
then the observed temperature of the breakout pulse is 

TbbmovI!,^^, yielding 

(13) 

If, in contrast, rj^o ^ 1, then the observed temperature 
is much higher than Tbb and is given by TBB,bo'r]bo^~'^ , 
where ^ accounts for the contribution of inverse comp- 
ton and its value is discussed below. Figure [T] depicts 
the observed breakout temperature as a function of the 
breakout velocity for tbo = 20 d. It is evident that Tsoft,bo 
depends strongly on deviation from thermal equilibrium, 
rising by more than an order of magnitude when Vbo is 
increased by a factor of 2. For example, for Vbo,9 = 0.3 
we find Tsoft,bo ~ 10^ K while if ff,o,9 = 1-5 then 
Tsoft,bo w 5 X 10^ K. The value of 77 that we calculate 
here0 fEauation [T2)) implies that the critical breakout ve- 

^ We assume here that the Thompson cross section is larger than 
the free- free absorption cross section and that r > 1. This is true 
while the soft component is dominant 

® Note that calculating r] for the configuration of a wind breakout 
is slightly different than in the steady state planar case discussed 
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locity, above which thermal equilibrium cannot be main- 
tained, is ^ 7000 kms^^. However, ry is calculated only 
to within an order of magnitude, and more importantly, 
the coefBcient in Equation [12] may be smaller by about 
an order of magnitude in high metallicity environment 
(see discussion in NSIO). Since the critical velocity for 
thermal equilibrium depends only weakly on that coef- 
ficient (to the power of —1/4), the critical velocity does 
not vary much between different systems. Nevertheless, 
since the velocity range that we consider here is exactly 
at the point where thermal equilibrium is marginal, we 
expect some of the breakouts in this range to be in equi- 
librium with Tsoft,bo ^ 10** K and some to be out of 
equilibrium with Tsoft,bo m the soft X-ray range. 

The evolution of Tbb and rjs (both measured at the 
shock radius) following the breakout is 



-BE. 



s{t) OC 



independent of n, and 
77,(t)oc^ocTP-5 



P 



BB 




(14) 



(15) 



The evolution of the soft component while thermal equi- 
librium is maintained is: 



6/11 




Tsoftit,ris < 1) cx Tbb,sV'. 



(16) 

Thus, while 77^ < 1, the observed temperature decreases 
slowly at t < tsp. It drops linearly with t during the 
snowplow phase for as long as ksTsoft ^ 1 eV, beyond 
which recombination takes place and the temperature 
drop stops. When rjs > 1 (deviation from thermal equi- 
librium) : 



if 1 
Tsoft{t,r]s > 1) OC Tbb.s-^ « 




(17) 



This equation is implicit since the factor ^ depends on 
Tsoft (see below) and its evolution in time is not specified 
explicitly. ^ is a logarithmic factor which accounts for the 
contribution of inverse compton in bringing the temper- 
ature closer to thermal equilibrium (^ always > 1). It 
is discussed in length in NSIO and its value is '--^ 1 at 
low temperatures (< 10^ K) and ~ 10 at high temper- 
atures (> 10^ K). In the scenario that we discuss here 
the calculation of ^ is slightly modified compared to the 
one described in NSIO, where the softest photon that is 



bv lWeaved ^973) and IKatz eFall pOia V The reason is that in a 
steady state shock photons that are diffusing from the downstream 
back to the shock, over the shock crossing time, are produced by 
a mass that is larger by the shock compression ratio, 7 in this 
case, than the mass in the shock transition layer. In the case 
of a wind breakout the mass behind the shock is comparable to 
the mass in the shock transition layer at the time of breakout. 
Therefore rjto is larger by a factor of Ri7 compared to the one that 
is calculated for the same shock velocity in steady state. As a result 
the radiation behind the wind breakout shock falls out of thermal 
equilibrium already when the shock exceeds 7000 kms~^ instead 
of 15,000 kms~^ in the steady state case. 



thermally coupled to the electrons is set by the require- 
ment that its energy is doubled before it is absorbed by 
free-free (see their Equation 11). Here, where the optical 
depth of the wind can become low enough, the limiting 
factor may become instead the low Compton y-parameter 
of the wind. In our calculations below we use both re- 
quirements (high enough y-parameter and no free-free 
absorption) to find the softest photons that are coupled 
by IC to the electrons. As r decreases the y-parameter 
becomes the limiting factor and the contribution of IC 
to coupling diminishes, so that Tsoft{t) — > TbbV^j a 
very rapid temperature rise before the snowplow phase. 

We note that in our calculation of T^oft we assume that 
the photons and the electrons have a single temperature 
at the region where Tsoft is determined. This is always 
true when 77 < 1, but it may be broken in some cases, 
especially when the shock velocity is low and r <C Tbo- If 
this happens then Tsoft is higher than the value that we 
calculate here. 

Finally, our analysis of Tsoft is not applicable any- 
more when the first of the following three events takes 
place: (i) recombination becomes important, (ii) the op- 
tical depth drops to the point where diffusion is not 
a good approximation (iii) the hard component carries 
most of the luminosity. 

3.2. Hard component 

The soft component is accompanied by a harder com- 
ponent, which is generated by the free-free emission 
from the hot electrons behind the shock. At breakout, 
the main cooling source of these electrons is IC of soft 
photons, yielding a free -free to IC emissivity ratio (c.f. 
iChevalier fc IrwinllMl : 



-IC.bo 



3/2 



A^/ksTee 



rad 



(18) 



where rig is the electron density, a is the fine structure 
constant, and trad ~ pv^ /2 at breakout. After the break- 
out, as long as T > 1, 

=.// 



1 



yi/2^2^ 



t t < tsp 
t>tsp 



(19) 



Therefore IC remains the dominant cooling source for a 
long time after the breakout, suppressing the emission 
of hard photons. But even the hard photons that are 
emitted at Tg are not seen directly by the observer as 
long as the wind optical depth is significant. Compton 
losses during their diffusion towards the observer limit 
their energy to meC^/r^ and farther reduces the lumi- 
nosity of the h ard component by a factor r'^Te/meC^ 
(jChevalier fc Irw in 2012). Therefore the observed tem- 
perature of the hard photons is 



7hard(i) 

and its luminosity is: 



r li\ T li\ -^hard (t) 

ihard(i) L{t) ^ mm 



Te{t) 



■ irff 



(20) 



(21) 



Equations [TOl [TSl [20l and Tbo = c/vbo imply that at break- 
out the hard component luminosity is always 

ihard.bo ~ lO^'^Lbo, (22) 
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and with Equations [71 and 1191 we find that at t < tgp, 
while -Lhaid ^ i it evolves as : 

ihard(i < isp) «i (23) 

If the hard component is still subdominant during the 
transition to the snowplow phase, or if the breakout takes 
place during this phase, then its fraction out of the to- 
tal luminosity rises quickly and multiplying this fraction 
with L{t): ' 

Li,,,d{t>tsp)(xt°~t^, (24) 

until it becomes comparable to the bolometric luminos- 
ity. 

3.3. Criterion for fast cooling 

We assumed above that the shocked plasma is cool- 
ing fast (within a dynamical time). Here we find when 
this assumption holds. A sufficient requirement for fast 
cooling is that free-free emission cools the hot electrons 
over a dynamical time. In case the compton y-parameter 
drops below unity, it is also a necessary condition. The 
free-free cooling time at breakout is much shorter than 
the dynamical time for our considered range of V{,o- 

t^f 1 

cool, bo ^ ^ 4 rp-1/2 ,r.r\ 

tbo ^0 

Past breakout, it evolves as 

^// f n = l2 

cx v'^n-^T-^/^-^ oc I n^7 (26) 

Thus, if the shock is cooling fast at tsp then it continues 
to cool fast also later. Requiring t^^l^^{tsp) /tsp < 1 and 
taking oc t'^ '^ {n = 12) all the way from breakout 

to tsp shows that the shock enters the snowplow phase 
cooling fast as long as t^o ^ Gvl^^g AI'^q^ E^-^^^ d. This 
calculation over estimates the threshold coefficient by a 
factor of ~ 3 since n varies from n < 12 at breakout to 
n w 4 at igp. Therefore, in the range of breakout veloci- 
ties that we consider, the shock is always fast cooling for 
breakout times that are days or longer. 

4. LIGHT CURVES OF THE DIFFERENT REGIMES 

The light curve and spectral evolution depend on three 
time scales. The first two are tbo and tsp, which affect 
the hydrodynamics. The third timescale, which we de- 
note thard, marks the time at which the hard component 
becomes the dominant one, thus affecting the spectral 
evolution. For breakout times of a day or longer, the or- 
der of appearance of these three characteristic timescales 
depends primarily on the breakout time tbo- 

Since at breakout the hard component is always a small 
fraction of the total luminosity, tbo < ^hard. Therefore, 
the evolution depends on the location of tsp with re- 
spect to tbo and thard- One regime is tsp < tbo < ^hard- 
The criterion for this regime was found in section [21 
When tbo < tsPj there are two additional regimes, sep- 
arated by the limiting case of tsp — ^hard- In order to 
find the separating criterion we use Equation [2TJ t^ard 
is the earliest time that satisfies Lhard ~ L, namely 
that both m mf,c^ /k^Tg and e-^^ « e^'~^ are satis- 
fied. Setting tsp = ihard the condition sa nieC^ /kpTe, 




t [days] 



Fig. 2. — A breakout within the snowplow regime, with = 70 
d and ng = 0.4 (matching £51 = 1.65 and m{,Q = IOMq), show- 
ing a very bright flash in the optical-UV, with a comparable rise 
time (not seen in the figure) and fall time. The bolometric lu- 
minosity and the temperature of the soft component {solid lines) 
decline as t~^-^ and t~^-^ respectively. Temperature decline stops 
around 6000 K, where recombination starts playing an important 
role (where the temperature line becomes dotted). The hard com- 
ponent luminosity and temperature {dashed lines) are strongly sup- 
pressed around the breakout (see text). > keV photons are ob- 
served only more than a year after the explosion. The breakout 
time and velocity are chosen to fit the observation of SN 2006gy. 

which is the last to be satisfied for typical parameters, 
implies r « 5EI(^ M-^^q^^^ , where we used v{tsp) ~ 
y^2E/AIej and Equation [TOl Using Equation [5l and set- 

^ /2 1/2 

ting T = 5Ec^[ A/]^Q one finds that tsp ~ ^hard when 
tbo,d ~ 20M^^'^''E^° '^^. We thus obtain the following 
three regimes with respect to tsp: 

tsp < tbo tbo^i > 80MiV-E5-i°-2^ 

tbo < tsP < fhard 20MiQ^ Er^^^'^^ < tbo.d ^ 80A/io^^ -B5 j^"'^ 

tbo < ihard < tsP 1 < tbo.d < 20Miq'''^ E ^-[''^^ 

(27) 

These regimes roughly match initial (breakout) n values 
of 4, 7 and 12, respectively. Below we discuss the lumi- 
nosity and spectral evolution in each of these regimes. 

4.1. Late breakout (tsp < tbo) 

tbo,d > 80AfiV5^5l"-25 

If the progenitor went through an extreme mass loss 
episode just prior to the explosion, the mass in the wind 
can be larger than the ejecta mass. In such a scenario the 
breakout may take place weeks to months after the explo- 
sion, near the time that, or even after, the reverse shock is 
exhausted. The entire explosion energy is released then 
around tbo- If the ejecta is several solar masses with 
~ 10^^ erg, the breakout velocity is 3000 — 5000 kms~^ 
and a very bright fiash, with a comparable rise time 
and fall time, is seen in the form of optical-UV photons. 
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with a typical temperature of ~ 1 — 5 x 10^ K. Follow- 
ing breakout, the luminosity and temperature of the soft 
component decline as t~^-^ and t~^'^ respectively. The 
temperature decline stops around 6000 K, where recom- 
bination starts playing an important role. The breakout 
hard component is strongly suppressed, (~ 10"'' of the 
total luminosity) by the IC cooling dominance and the 
high wind optical depth. Moreover, its temperature is 
limited during the breakout to < 0.5 keV by Compton 
loses during the diffusion through the unshocked wind 
([chevalier & Irwin 2012). The hardness and luminosity 
of the hard component rise following the breakout, but a 
significant emission above 1 keV starts only at ~ lOtto- 
Figure [2] depicts the luminosity and temperature evolu- 
tion of a breakout near the snowplow phase with tto = 70 
d and VQ^bo = 0.4 where the evolution discussed above is 
seen. The breakout time scale and velocity fit t he ob- 
servations of SN 2006gy (lOfek et al. 2007:' Smith et al.l 
[200l . which was suggested bv lChevalier fc IrwinI ()2011f ) 
to be a shock breakout through a wind. The luminos- 
ity and temperature of the dominating soft component 
of SN 2006gy are recovered with these ti,o and Vbo- Sim- 
ilarly to the conclusion of Chevalier & Irwin (2OI4), we 
find a strongly suppressed X-ray luminosity, with keV 
photons appearing only more than a year after the ex- 
plosion (assuming that the wind is not ended at an earlier 
time). 

Luminous hard X-ray emission is expected only if the 
breakout velocity is high, which requires either a low 
ejecta mass ~ Mq or a very high energy explosion ^ 10^^ 
erg. When the breakout velocity is high and Uo is low 
(< 30) the hard component releases most of its energy 
at > 1 keV photons soon after the breakout. The lumi- 
nosity of SNe in this phase is so high that X-rays may be 
bright enough for detection even though they carry only 
a small fraction of the total energy. 

Note that already during the breakout the mass of the 
wind is comparable to the mass of the ejecta, and the 
mass accumulated by the shock grows as t^/^. Therefore, 
if the wind mass is not much larger than the ejecta mass, 
the shock encounters the end of the massive wind not too 
long after tbo- 

4.2. Early breakout, (tto < thard < tsp) 

1 < hoM < 20MiV5iJ-"-25 

If breakout occurs early the shock velocity is larger, rang- 
ing from 5,000-10,000 kms^^ for typical SNe. The break- 
out emission is therefore at a marginal thermal equilib- 
rium and a breakout temperature in the range 10"*^— 10^ K 
is expected for SNe with typical energy and ejecta mass. 
Figure [3] depicts the luminosity and temperature of the 
soft and hard components in a breakout with tf,o = 10 
d and Vbo = 5,000 kms~^. The soft component is in 
thermal equilibrium during breakout, rjto — 0.3, result- 
ing in Tsoft,bo ~ 25, 000 K. Following the breakout 77s 
increases roughly linearly with time, implying that it 
becomes harder for the radiation to maintain thermal 
equilibrium. As a result Tsoft declines very slowly until 
day 25, when thermal equilibrium cannot be maintained 
anymore, and rises quickly afterwards. The luminosity 
decreases slowly during this time. The flux in a fixed 
optical/UV band is expected to drop slowly at first and 
then much faster. Such a breakout may explain SN PTF 




10' 10' 10' 

t [days] 



Fig. 3. — Bolometric luminosity {solid line), soft component tem- 
perature (solid line), and luminosity and temperature of the hard 
component (dashed lines) from an early breakout of t^o = 10 d and 
^bo = 5i 000 kms~^. tgp is calculated for E51 = 1 and Mio = 1.25. 
The temporal decay of the bolometric luminosity is shifted from 
(n=12) to t-°-^ (n=7) at O.ltsp, to account for the flattening 
of the ejecta profile once m approaches M^j . The soft component 
is in thermal equilibrium during breakout, rji,^ = 0.3, resulting in 
Tsoft,ho ~ 25, 000 K. Tgoft decreases slowly before thermal equi- 
librium is broken, at day 25, and rises rapidly afterwards (the soft 
temperature line is cut at t = 3, where diffusion approximation 
does not hold an ymore). S uch a breakout may explain SNe like 
PTF 09uj (Ofek et al.1120101 '1. The hard component, very faint and 
soft at first, becomes brighter and harder quickly and X-rays/soft 
gamma-rays dominate the flux at day 200. Such events make ex- 
cellent candidates for X-ray searches ~100-300 d after the SN ex- 
plosion. 

09uj, as suggested bv lOfek et all (|2010l ). The breakout 
temperature and luminosity of the case depicted in figure 
[3] are consistent with the observations (jOfek et al.l HolOl. 
and so does the general trend of the fading i?-band fiux, 
which becomes faster with time. The hard component is 
very faint at first (only ~ lO"'^ of the bolometric lumi- 
nosity) and its temperature is limited to ^ 0.5 keV. But 
it becomes brighter and harder quickly, with hard X-rays 
and soft gamma-rays dominating the fiux at day 200. 

In the general case the breakout thermal coupling pa- 
rameter Tyho may be smaller or larger than the one in 
the case depicted in figure |3l If it is smaller then ther- 
mal equilibrium is kept for longer time, while if it is 
larger then the breakout temperature can be significantly 
larger and it is starts to climb fast immediately after 
breakout. The hard component, as always, starts with 
^ 10""* of the total luminosity and rises quickly to dom- 
inate the luminosity at ^ 20tf,o- Since the luminosity 
decays slowly during this time, these events are good 
candidates for X-ray searches ^100-300 day after the SN 
explosion, when X-rays dominate the luminosity and be- 
fore the snowplow phase starts. Finally, the luminosity 
increases with breakout time. Therefore, breakouts with 
larger tbo (which are still in this regime) produce brighter 
X-ray events. 

4.3. Intermediate breakout time, (tbo < tsp < ^hardj 
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t [days] 



Fig. 4. — An intermediate regime breakout, of t(,o = 30 d and 
^bo = 5,000 kms~^. tgp is calculated for Esi = 1 and Afio = 1. 
The transition to the snowplow phase happens at day 230, while 
the hard component gains dominance only after this transition, 
reaching a lower luminosity and softer X-ray temperatures com- 
pared to those of the early regime. If the breakout is early enough, 
so that thard ~ ^SPi such events produce the brightest X-ray lu- 
minosity among all the cases. The soft component temperature 
decline stops around 6000 K due to recombination (where the tem- 
perature line becomes dotted). Notations are similar to figures [2] 
and|3] 

20MiV5£;-o.25 < ^^^^ < 80MiV'£^5"i° '' 

Figure |4] presents the luminosity and temperature of an 
intermediate case, where the breakout takes place before 
the transition to the snowplow phase, but the hard com- 
ponent becomes dominant only after this transition. The 
properties of the breakout and the emission that follows 
soon after it are similar to those described in the early 
breakout case. The difference starts at the transition to 
the snowplow phase, where the temperatures of both the 
soft and the hard components change direction and start 
to decrease with time. The main difference from the early 
breakout case is that the hard component does not dom- 
inate the luminosity at tsp and by the time that it does, 
the luminosity is already dropping fast. If the breakout is 
early enough and ihard ~ ^sp, then the X-ray luminosity 
is high. In fact these events produce the brightest X-ray 
luminosity among all the cases discussed in this paper. 
If, however, thard ^ tsp, then the bolometric luminosity 
can drop significantly before X-rays dominate, resulting 
in a much dimmer signal. 

5. SUMMARY 

We examine the emission from a shock breakout 
through a thick wind. We restrict our exploration to 
breakout time scales that are days or longer, which are 
easier to detect, and velocities in the range 3,000-15,000 
kms^^, as expected for typical SNe. We also consider 
only a standard wind profile pw esc r~^, that is extended 
beyond the breakout radius. Following the breakout, 
the internal energy in these systems is dominated by 
the wind material that is heat ed by the forwar d shock 
(fChevalier fc Irwin ,2011.: .Balberg fc Loebll2011[) . which 



makes a transition from a radiation mediated shock to a 
collisionless shock at the breakout ()Katz et al.ll201l"bl) . 
Our main conclusions are listed below. 

(i) The breakout and following emission are very 
luminous. In the regime that we explore, the shocked 
plasma is always in the fast cooling regime (see the 
exact condition in section [3.31) . converting all the shock 
luminosity into radiation. Therefore, if the wind mass 
is comparable to, or larger than, the ejecta mass, then 
all the SN energy is radiated away. If the wind mass, 
Mw is smaller than Mej then the energy radiated during 
the interaction is ^ EM^/Mej. This re sult supports 
the suggestion by iChevalier fc IrwinI (|2011[ ) that at least 
some ultra-luminous type Iln SNe are such breakouts. 

(ii) The post-breakout emission is composed of two 
spectral components - soft and hard. The electrons 
behind the collisionless shock are heated to temper- 
atures of - min [60, 200t;^] keV (jKatz et all l2011bD . 
The soft component is radiated by the IC emission of 
these electrons over photons produced by the heated 
unshocked wind, while the hard component is generated 
by the free-free emission of the shocked electrons. 

(iii) The soft component is in marginal thermal equi- 
librium during the breakout. Thus, the breakout 
temperature can be in the range ^ 10'' — 10^ K, where 
slower breakouts, w 3,000 kms~^ will be in the lower 
side of this range while faster breakouts w 15, 000 
kms~^ will populate the high end. In the latter case, 
the breakout X-ray emission of the soft component 
exponential tail may very well dominate over the hard 
component, and may possibly be detectable. 

(iv) The hard component is suppressed as long as IC 
is more efficient than free-free and the wind optical 
depth is large (Chevalier & Irwin 2012). We find that 
during breakout the hard component luminosity is 
always a fraction of ~ 10~^ of the soft one. This result 
is independent of breakout time and velocity (it is also 
independent of the wind density profile as long as it is 
not terminated at the breakout radius). The fraction 
of the hard component out of the total luminosity 
rises quickly after breakout. It becomes dominant at 
- 10- 50iho. 

The post-breakout evolution of the two components 
can be divided into three regimes according to the break- 
out time. The behavior in each of the regimes is as fol- 
lows: 

Late breakout (typically tbo > 80 d): Breakouts in the 
snowplow regime (i.e., when the shocked wind mass is 
already larger than M^j) provide the brightest emission 
among all the regimes. All the explosion energy is re- 
leased as optical-UV photons and no X-ray emission is 
expected during the breakout. X-rays are seen only long 
after the breakout in the case of regular SNe or soon af- 
ter the breakout in the case of super-energetic SNe (high 
velocity breakout). At the peak of the hard component, 
around 20tbo, X-rays carry only ^ 10"^ of the breakout 
energy, in the form of ^ 1—10 keV photons. Follow- 
ing the breakout, the temperature of the soft component 
drops roughly as t~^'^ until recombination starts. The 
bolometric luminosity drops as (also after recombi- 
nation takes place), at first dominated by the soft com- 
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ponent, which drops quickly once the hard component 
becomes dominant. 

Early breakout (typically 1 < ti,o < 20 d): The break- 
out radiation is in marginal thermal equilibrium (as dis- 
cussed above), and following the breakout it becomes 
harder to maintain thermal equilibrium. As a result, if 
the soft component is in thermal equilibrium at breakout, 
it keeps roughly a constant temperature until it falls out 
of thermal equilibrium, and then the temperature climbs 
quickly. If the soft component is out of thermal equilib- 
rium already at breakout, then its temperature starts to 
rise immediately. The hard component becomes domi- 
nant before the transition to the snowplow phase, while 
the velocity of the shock is high and the energy behind 
the shock is still rising (the luminosity at this point is 
comparable to that of the breakout). As a result, these 
bursts are predicted to produce a bright hard X-rays to 
soft gamma-rays emission on time scales of ^ 20tbo- 
Intermediate breakout time (typically 20 d < tbo < 80 d) : 
The breakout takes place before the transition to the 
snowplow phase, but the hard component becomes dom- 
inant only after this transition. The emission following 
the breakout is similar to the early breakout case until 
the transition to the snowplow phase, which takes place 
before the hard component dominates. As a result, the 
peak X-ray luminosity is lower than the breakout lumi- 
nosity. 

In our calculations we ignore the radius at which the 
thick wind ends. Once the wind density drops abruptly 
the forward sh ock becomes inefficient a nd the luminosity 
fades quickly (jChevalier fc Irwinll2011|) . Since the mass 
of the progenitor at birth is limited, and the mass in the 
wind needed for an extended breakout is considerable, it 
is expected that the light curves we present will be termi- 
nated at some point. For example, in order for the light 
curve to remain bright until tsp the wind mass must be 
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> Mej- Even when > M^j the collected mass is 
oc i^/^ after tsp, limiting the lifetime of the bright emis- 
sion. In terms of the prospects of observing X-rays, this 
consideration boosts the attractiveness of earlier break- 
outs, since the mass that early breakout events collect by 
the time that the hard component dominates is lower. 

We do not address the extinction of soft X-rays due 
to photoabsorption by pa rtially ionized unshocked wind 
(Chevalier & Ir^ [20121 ). This process may suppress 
softer-end photons of the hard component as well as hard 
photons of the soft component, if the latter is far from 
thermal equilibrium. If photoabsorption is important 
then the observed flux of soft X-rays is lower than the 
one that we predict. 

Finally, our treatment ignores the contribution of emis- 
sion from ejecta mass before it is shocked by the reverse 
shock. Once the Compton y-parameter of the wind drops 
below unity, this emission becomes similar to the emis- 
sion from a typical SN without a thick wind (e.g., typical 
IIP, lb and Ic SNe) . The emission from the forward shock 
that we consider here always outshines this emission (in 
terms of bolometric luminosity). But, at late times, when 
most of the forward shock emission is in hard photons, 
the contribution from inner layers may dominate the IR- 
optical emission. 

To conclude, breakouts through thick winds produce 
very bright SNe with a potential for detection across the 
entire electromagnetic spectrum. The brightest phase 
is the breakout, which is always dominated by optical- 
UV emission with a negligible X-ray emission. The X- 
rays and even soft gamma-rays luminosity rises after the 
breakout, and reaches its peak after ^ 10 — 50 ^bo- 
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